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Abstract 



With a large integrated luminosity expected at the Tevatron, a next-to- 
leading order (NLO) calculation is no longer sufficient to describe the data 
which yield the precision measurement of etc. Thus, we extend the 

CoUins-Soper-Sterman resummation formalism, for on-shell vector boson pro- 
duction, to correctly include the effects of the polarization and the width of 
the vector boson to the distributions of the decay leptons. We show how to 
test the rich dynamics of the QCD multiple soft gluon radiation, for example, 
by measuring the ratio Rcss = — ^ — . {Qt is the transverse momentum 
of the vector boson.) We conclude that both the total rates and the distri- 
butions of the lepton charge asymmetry predicted by the resummed and the 
NLO calculations are different when kinematic cuts are applied. 
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I. INTRODUCTION 



Quantum Chromodynamics (QCD) is a field theory that is expected to explain all the ex- 
perimental data involving strong interactions either perturbatively or non-perturbatively 
Consider the weak boson {W^ and Z^) production at a hadron collider, such as the Teva- 
tron. In the framework of QCD the production rate of the weak bosons is calculated by 
multiplying the constituent cross section (the short- distance or the perturbative physics) by 
the parton luminosities (the long-distance or the non-perturbative physics) [Q. This pre- 
scription of theoretical calculation was proven to the accuracy of (^(l/Q^) and is known as 
the factorization theorem of QCD |^. (Q is the mass of the vector boson.) Since we do not 
yet know how to solve QCD exactly, we have to rely on the factorization theorem to sepa- 
rate the perturbative part from the non-perturbative part of the formalism for any physical 
observable. The short distance contribution can be calculated perturbatively order by order 
in the strong coupling as- The long distance part has to be parametrized and fitted to the 
existing data so that it can later be used to predict the results of new experiments. There- 
fore, theoretical predictions that are compared to experimental data always has to invoke 
some approximation in the calculations based upon QCD. We refer to different prescriptions 
of calculations to be different models of theory calculations which all originate from the one 
and only QCD theory. For instance, to improve theory predictions on the event shape, such 
as the transverse momentum (Qt) distribution of the weak boson, the commonly used theory 
model is the event generator, e.g., ISA JET |], PYTHIA §, or HERWIG 0. The event 
generator can also provide information on the particle multiplicities or the number of jets, 
etc. 

However, as discussed above, different models of calculation make different approxima- 
tions. Hence, a model can give more reliable theory predictions than the others on some 
observables, but may do worse for the other observables. A few more examples are in order. 
To calculate the total production rate of a weak boson at hadron colliders, it is better to use 
a fixed order perturbation calculation, and a higher order calculation is usually found to be 
more reliable than a lower order calculation because it is usually less sensitive to the choice 
of the scale for calculating the parton distribution functions (PDF) or the constituent cross 
section (including the strong coupling constant as). The former scale is the factorization 
scale and the latter is the renormalization scale of the process. Unfortunately, a fixed or- 
der perturbation calculation cannot give reliable prediction of the distribution of Qt when 
Qt is small. On the contrary, an event generator can give more reliable prediction for Qt 
distribution in the small Qt region, but it usually does not predict an accurate event rate. 
The general feature of the above two theory models is that a fixed order calculation is more 
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reliable for calculating the event rate but not the event shape, and an event generator is 
good for predicting the event shape but less reliable for the event rate. 

In this paper, we discuss another model of theory calculation that can give reliable predic- 
tions on both the event rate and the shape of the distributions. Specifically, we are interested 
in the distributions of the weak bosons and their decay products. This model of calculation 
is to resum a series of large perturbative contributions due to soft gluon emission predicted 
by the QCD theory. We present the QCD resummation formalism for calculating the fully 
differential cross section of the hadronically produced lepton pairs through electroweak (EW) 
vector boson production and decay: hih2 * (ii(i2)X. We focus our attention on the 

Tevatron though our calculation is general and applicable for any hadronic initial state hih2 
and any colorless vector boson. For instance, the vector boson V can be one of the standard 
model (SM) electroweak gauge bosons or the virtual photon 7* (for producing the 
Drell-Yan pair), or some exotic vector boson such as Z' and W in the extended unified gauge 
theories. 

At the Tevatron, about ninety percent of the production cross section of the and 
bosons (with mass Q) is in the small transverse momentum {Qt) region, where Qt ^ 20 
GeV (hence <^ Q"^)- In this region the higher order perturbative corrections, dominated 
by soft and coUinear gluon radiation, of the form Sn^i I^m=o^ ""^m «sln™'(QT/Q^); 
substantial because of the logarithmic enhancement [0. („fm are the coefficient functions 
for a given n and m.) These corrections are divergent in the Qt — > limit at any fixed 
order of the perturbation theory. After applying the renormalization group analysis, these 
singular contributions in the low Qt region can be resummed to derive a finite prediction for 
the Qt distribution to compare with experimental data. It was proven by Collins and Soper 
in Ref. p that not only the leading logs |l9|JT0(] but all the large logs, including the sub-logs 
in the perturbative, order-by-order calculations can be resummed for the energy correlation 
in e^e~ collisions. 

For the production of vector bosons in hadron collisions two different formalisms were 
presented in the literature to resum the large contributions due to multiple soft gluon radia- 
tion: by Altarelli, Ellis, Greco, Martinelli (AEGM) [|ll| ; and by Collins, Soper and Sterman 



(CSS) 0. The detailed differences between these two formalisms were discussed in Ref. |]T2 



It was shown that the AEGM and the CSS formalisms are equivalent up to the few highest 
power of ln(Qy/(5^) at every order in as for terms proportional to Qt^, provided as in the 
AEGM formalism is evaluated at h'^/h'^ rather than at Q^. A more noticeable difference, 
except the additional contributions of order Q~^ included in the AEGM formula, is caused 
by different ways of parametrizing the non-perturbative contribution in the low Qt regime. 
Since the CSS formalism was proven to sum over not just the leading logs but also all the 
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sub-logs, and the piece including the Sudakov factor was shown to be renormalization group 
invariant 0, we only discuss the results of CSS formalism in the rest of this paper. 

With the increasing accuracy of the experimental data on the properties of and Z° 
bosons at the Tevatron, it is no longer sufficient to only consider the effects of multiple 
soft gluon radiation for an on-shell vector boson and ignore the effects coming from the 
decay width and the polarization of the massive vector boson to the distributions of the 
decay leptons. Hence, it is desirable to have an equivalent resummation formalism |]TE[ for 
calculating the distributions of the decay leptons. This formalism should correctly include 
the off-shellness of the vector boson (i.e. the effect of the width ) and the polarization 
information of the produced vector boson which determines the angular distributions of the 
decay leptons. 

In the next section, we give our analytical results for such a formalism that correctly 
takes into account the effects of the multiple soft gluon radiation on the distributions of 
the decay leptons from the vector boson. In Section III, we discuss the phenomenology 
predicted by this resummation formalism. To illustrate the effects of multiple soft gluon 
radiation, we also give results predicted from a next-to-leading order (NLO) calculation. 
As expected, the observables that are directly related to the transverse momentum of the 
vector boson will show large differences between the resummed and the NLO predictions. 
These observables are the transverse momentum of the leptons from vector boson decay, the 
back-to-back correlations of the leptons from Z° decay, etc. The observables that are not 
directly related to the transverse momentum of the vector boson can also show noticeable 
differences between the resummed and the NLO calculations if the kinematic cuts applied 
to select the signal events are strongly correlated to the transverse momentum of the vector 
boson. Section IV contains our detailed discussion. 

Since this Qt resummation formalism only holds in the Collins-Soper (CS) frame [l^, we 
give the detailed form of the transformation between a four-momentum in the CS frame (a 
special rest frame of the vector boson) and that in the laboratory frame (the center-of-mass 
frame of the hadrons hi and in Appendix A. In Appendix B the analytical expression 
for the NLO results, of 0{as), are given in D = 4 — 2e dimensions. Appendix C contains 
the expansion of the resummation formula up to 0{as)- Appendix D lists the values of A, 
B, and C functions (cf. Sec. II) used for our numerical calculations. 

We note that the resummation formalism presented in this paper can be applied to any 
processes of the type hih2 —>■ V{—>- ^1^2)-^? where ^ is a color neutral vector boson which cou- 
ples to quarks and leptons via vector or axial vector currents, that is = ej, W^, Z^, W, Z', 
etc. Throughout this paper, we take V to be either or Z^ bosons, unless specified oth- 
erwise. 
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II. THE RESUMMATION FORMALISM 



To derive the resummation formalism, we use the dimensional regularization scheme to 
regulate the infrared divergencies, and adopt the canonical-75 prescription to calculate the 
anti-symmetric part of the matrix element in D-dimensional space-time.^ The infrared- 
anomalous contribution arising from using the canonical-75 prescription was carefully han- 



dled by applying the procedures outlined in Ref. |T8[ for calculating both the virtual and the 
real diagrams .0 

The kinematics of the vector boson V (real or virtual) can be expressed in terms of 
its mass Q, rapidity y, transverse momentum Qt, and azimuthal angle (pv, measured in 
the laboratory frame (the center-of-mass frame of hadrons hi and h2). The kinematics of 
the lepton ii is described by 6 and 0, the polar and the azimuthal angles, defined in the 
Collins-Soper frame which is a special rest frame of the ^-boson [|l9l. (A more detailed 



discussion of the kinematics can be found in Appendix A.) The fully differential cross section 
for the production and decay of the vector boson is given by the resummation formula in 
Ref. [H: 



f da{hih2-^V{^ hl2)X) \ ^ _^ 

[dQ^ dydQld(f)vd cos ed(j)) GGn^S {Q"^ - M^f + Q^VH 



X 



(27r) 



j d^e'"^^-' J2 Wfkib*, Q, xi,X2, 9, 0, Ci, C2, C3) Wjl'^ib, Q, xi, X2) 



+ YiQT,Q,xi,X2,9,(P,C^)}. (1) 

In the above equation the parton momentum fractions are defined as Xi = e^QI\fS and 
X2 = e~'^Q/\/S, where ^/S is the center-of-mass (CM) energy of the hadrons hi and h2. 
For V = or we adopt the LEP line-shape prescription of the resonance behavior. 
The renormalization group invariant quantity Wjj^{b), which sums to all orders in as all the 
singular terms that behave as Q^'^ x [1 or In (Q^/Q^)] for Qt — > 0, is 

Wf,{b, Q, xi, X2, e, 0, Ci, C2, Cs) = exp {-S{b, Q, Ci, C2)} \ V,k \' 

X { [(Cj„ ® fa/h^) (Xi) (Cfcfe ® fb/h^) (Xa) + ® fa/h,) (Xi) (^Cjb ® fb/h^) {X2) 



^In this prescription, 75 anti-commutes with other 7's in the first four dimensions and commutes 
in the others |16,17|. 



^In Ref. [^], the authors calculated the anti-symmetric structure function for deep-inelastic 
scattering. 
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+ 

X 



{Cja ® fa/ hi) (Xl) [Ckb ® fb/h2) ix2) - (C^^ ® fa/h^) (Xi) (Cjb ® fb/hi) (3^2 ) 

i9l-9l)ifl-fR)i2cose)}, (2) 



where ® denotes the convolution 



(3) 



and the Vjk coefficients are given by 



_ \ Cabibbo — Kobayashi — Maskawa matrix elements for V = ^ . 

''~\Sjk for\/ = Z0,7*- 

In the above expressions j represents quark flavors and k stands for anti-quark flavors. The 
indices a and b are meant to sum over quarks and anti-quarks or gluons. Summation on 
these double indices is implied. In Eq. we define the couplings fL,R and gL,R through the 
£1^2 and the qq'V vertices, which are written respectively, as 

«7m[/l(1 -75) + /r(1 + 75)] and - 75) + + 75)] . 

For example, for V = W~^,q = u, q' = d, ii = i/g, and £2 = e"*", the couplings 9^ = fl = 
GpM^/\/2 and 9r = fu = 0, where Gp is the Fermi constant. The detailed information 
on the values of the parameters used in Eqs. (|l]) and (|^) is given in Table |. The Sudakov 
exponent 5(6, Q, Ci, C2) in Eq. (|^) is defined as 



S{b,Q,C,,C2)= f 
Jc 



2/n2' 

(5) 



A {asm), Ci) In f ^ 1 + B [asifi), C,, 



The explicit forms of the A, B and C functions and the renormalization constants Ci (i=l,2,3) 
are summarized in Appendix D. 

In Eq. d^) the magnitude of the impact parameter h is integrated from to 00. However, 
in the region where h ^ I/Kqcd, the Sudakov exponent 5(6, Q, Ci, C2) diverges as the 
result of the Landau pole of the QCD coupling ^^(/i) at yU = Aqcd, and the perturbative 
calculation is no longer reliable. As discussed in the previous section, in this region of the 
impact parameter space (i.e. large b), a prescription for parametrizing the non-perturbative 
physics in the low Qt region is necessary. Following the idea of Collins and Soper the 
renormalization group invariant quantity Wji{b) is written as 



Wf,{b) = P^,s(6.)W^r (6) 
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V 


My (GeV) 


Tv (GeV) 


ah 


9R 


7 


0.00 
80.36 
91.19 


0.00 
2.07 
2.49 


gQfSu,/2 
5/(2^2) 
g{n-Qfsl)/{2c^) 


gQfSw/2 





TABLE I. Vector boson parameters and couplings to fermions. The f f'V vertex is defined as 
z7^[5fL(l— 75) + (7ij(l+75)] and = smOw {cw = cos Ow) is the sine (cosine) of the weak mixing angle: 
s\ii^{ey,{Mzo))j^ = 0.2315. Qf is the fermion charge (Q„ = 2/3, = -1/3, Q^- = -^,Qu = 0), 
and T3 is the eigenvalue of the third component of the SU{2)l generator (T3* = 1/2, = —1/2, 

= 1/2, Tf = -1/2). 

Here Wji^{b^:) is the perturbative part of Wji^{b) and can be reliably calculated by perturbative 
expansions, while W^^{b) is the non-perturbative part of Wjj.{b) that cannot be calculated 
by perturbative methods and has to be determined from experimental data. To test this 
assumption, one should verify that there exists a universal functional form for this non- 
perturbative function W^^{b). This is similar to the general expectation that there exists a 
universal set of parton distribution functions (PDF's) that can be used in any perturbative 
QCD calculation to compare it with experimental data. In the perturbative part of Wjj^{b), 

+ {b/b^accy ' 

and the non-perturbative function was parametrized by (cf. Ref. 0) 



Wjk (b, Qo, xi,X2) = exp 



-Fi{b) In ( ^ ) - Fj /h.ixi, b) - ^^//^^(xa, b) 



0. 



(6) 



where Fi, Fj/h^ and -Ffc/hj have to be first determined using some sets of data, and later can 
be used to predict the other sets of data to test the dynamics of multiple gluon radiation 
predicted by this model of the QCD theory calculation. As noted in Ref. 0, Fi does not 
depend on the momentum fraction variables xi or X2, while -Fj/hi and Fj^j^^ in general depend 
on those kinematic variables.^] The hi{Q^ /Q"^) dependence associated with the Fi function 
was predicted by the renormalization group analysis 0. Furthermore, Fi was shown to 
be universal, and its leading behavior (~ 6^) can be described by renormalon physics ||20|| . 
Various sets of fits to these non-perturbative functions can be found in Refs. pl|] and p2|. 



^Here, and throughout this work, the flavor dependence of the non-perturbative functions is ig- 
nored, as it is postulated in Ref. [0]. 
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In our numerical results in the next section, we use the Ladinsky-Yuan parametrization 
of the non-perturbative function (cf. Ref. p2|): 



Wj^ {b, Q, Qo, xi, X2) = exp 



-gib"^ - fl'2&^ln - fi-ifl-sfe In (100x1X2) 



(7) 



where gi = 0.1ll[];[J^ GeV^ ^2 = 0.58lJ]i GeV^ c/3 = -l-^-Vi GeV"^ and Qo = 1.6 GeV. 
(The value b^ax = 0.5 GeV^^ was used in determining the above gi's and in our numerical 
results.) These values were fit for CTEQ2M PDF with the canonical choice of the renormal- 
ization constants, i.e. Ci = C3 = 2e~'^'^ {'Je is the Euler constant) and C2 = 1. In principle, 
for a calculation using a more update PDF, these non-perturbative parameters should be 
refit using a data set that should also include the recent high statistics Z° data from the 
Tevatron. This is however beyond the scope of this paper. 

In Eq. (|ll), Wj^ sums over the soft gluon contributions that grow as or 
ln(Q^/Q^)] to all orders in as- Contributions less singular than those included in Wj% 
should be calculated order- by-order in as and included in the Y term, introduced in Eq. (||). 
This would, in principle, extend the applicability of the CSS resummation formalism to all 
values of Qt- However, as to be shown below, since the A, B, C, and Y functions are only 
calculated to some finite order in as, the CSS resummed formula as described above will 
cease to be adequate for describing data when the value of Qt is in the vicinity of Q. Hence, 
in practice, one has to switch from the resummed prediction to the fixed order perturbative 
calculation as Qt ^ Q- The Y term, which is defined as the difference between the fixed 
order perturbative contribution and those obtained by expanding the perturbative part of 
Wj^ to the same order, is given by 



y{Qt, q, xi, X2, e, <p, c) = f'^ f'^f^ 

Jxi ^1 Jx2 C,2 n=l 



asiC^Q) 



TT 



xfa/h, (^1, C,Q) R'2\Qt, Q, ^, 0) /,/h.(6, C4Q), (8) 

where the functions i?^^^ contain contributions less singular than Q^"^ x [1 or In (Q^/Q^)] 
as Qt 0. Their explicit expressions and the choice of the scale C4 are summarized in 
Appendix D. 

Before closing this section, we note that the results of the usual next-to-leading order 
(NLO), up to 0{as), calculation can be obtained by expanding the above CSS resummation 
formula to the as order, which includes both the singular piece and the Y term. Details are 
given in Appendices B and D, respectively. 
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III. PHENOMENOLOGY 



As discussed above, due to the increasing precision of the experimental data at hadron 
coUiders, it is necessary to improve the theoretical prediction of the QCD theory by including 
the effects of the multiple soft gluon emission to all orders in as- To justify the importance 
of such an improved QCD calculation, we compare various distributions predicted by the 
resummed and the NLO calculations. For this purpose we categorize measurables into two 
groups. We call an observable to be directly sensitive to the soft gluon resummation effect 
if it is sensitive to the transverse momentum of the vector boson. The best example of such 
observable is the transverse momentum distribution of the vector boson {da /dQT)- Likewise, 
the transverse momentum distribution of the decay lepton (da/dp^) is also directly sensitive 
to resummation effects. The other examples are the azimuthal angle correlation of the two 
decay leptons {A(j)^'^^^), the balance in the transverse momentum of the two decay leptons 
{Pt ~Pt)i ^he correlation parameter z = —pr^^ ■j?/^/[max(p'/^,p^^)]^. These distributions 
typically show large differences between the NLO and the resummed calculations. The 
differences are the most dramatic near the boundary of the kinematic phase space, such 
as the Qt distribution in the low Qt region and the A^^^^^ distribution near vr. Another 
group of observables is formed by those which are indirectly sensitive to the resummation 
of the multiple soft gluon radiation. The predicted distributions for these observables are 
usually the same in either the resummed or the NLO calculations, provided that the Qt is 
fully integrated out in both cases. Examples of indirectly sensitive quantities are the total 
cross section cr, the mass Q, the rapidity y, and Xp {= 2q^/\/S^ of the vector bosonQ, and 
the rapidity of the decay lepton. However, in practice, to extract signal events from 
the experimental data some kinematic cuts have to be imposed to suppress the background 
events. It is important to note that imposing the necessary kinematic cuts usually truncate 
the range of the Qt integration, and causes different predictions from the resummed and the 
NLO calculations. We demonstrate such an effect in the distributions of the lepton charge 
asymmetry A{y^) predicted by the resummed and the NLO calculations. We show that 
they are the same as long as there are no kinematic cuts imposed, and different when some 
kinematic cuts are included. They differ the most in the large rapidity region which is near 
the boundary of the phase space. 

To systematically analyze the differences between the results of the NLO and the re- 
summed calculations we implemented the 0{ag) (LO), the 0{as) (NLO), and the resummed 



^Here is the longitudinal-component of the vector boson momentum q^^ [cf. Eq. ( Al)] 
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FIG. 1. Total production cross section as a function of the parameter Q^"^^ (solid curve). 
The long dashed curve is the part of the 0{as) cross section integrated from Q^'^^ to the kinematical 
boundary, and the short dashed curve is the integral from Qt = to Q^^^ at 0{as)- The total 
cross section is constant within 10~^ % through more than two order of magnitude of Q^'^^. 

calculations in a unified Monte Carlo package: ResBos (the acronym stands for Resnmmed 
Vector i^oson production). The code calculates distributions for the hadronic production 
and decay of a vector bosons via hih2 —>■ V{~>- iii2)X, where hi is a proton and /i2 can be a 
proton, anti-proton, neutron, an arbitrary nucleus or a pion. Presently, V can be a virtual 
photon 7* (for Drell-Yan production), or Z^. The effects of the initial state soft gluon 
radiation are included using the QCD soft gluon resummation formula, given in Eq. (|l]). This 
code also correctly takes into account the effects of the polarization and the decay width of 
the massive vector boson. 

It is important to distinguish ResBos from the parton shower Monte Carlo programs like 
ISA JET II, PYTHIA [|, HERWIG 0, etc., which use the backward radiation technique || 
to simulate the physics of the initial state soft gluon radiation. They are frequently shown 
to describe reasonably well the shape of the vector boson distribution. On the other hand, 
these codes do not have the full resummation formula implemented and include only the 
leading logs and some of the sub-logs of the Sudakov factor. The finite part of the higher 
order virtual corrections which leads to the Wilson coefficient (C) functions is missing from 
these event generators. ResBos contains not only the physics from the multiple soft gluon 
emission, but also the higher order matrix elements for the production and the decay of the 
vector boson with large Qt, so that it can correctly predict both the event rates and the 
distributions of the decay leptons. 

In a NLO Monte Carlo calculation, it is ambiguous to treat the singularity of the vector 
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boson transverse momentum distribution near Qt = 0. There are different ways to deal with 
this singularity. Usually one separates the singular region of the phase space from the rest 
(which is calculated numerically) and handles it analytically. We choose to divide the Qt 
phase space with a separation scale . We treat the Qt singular parts of the real emission 
and the virtual correction diagrams analytically, and integrate the sum of their contributions 
up to ^. If Qt < Qt"^ we assign a weight to the event based on the above integrated 
result and put it into the Qt = bin. If Qt > Qt'^, the event weight is given by the usual 
NLO calculation. The above procedure not only ensures a stable numerical result but also 
agrees well with the logic of the resummation calculation. In Fig. |l| we demonstrate that the 
total cross section, as expected, is independent of the separation scale Qr^^ in a wide range. 
As explained above, in the Qt < Qt^^ region we approximate the Qt of the vector boson 
to be zero. For this reason, we choose Q^^ as small as possible. We use Qt^^ = 0.1 GeV 
in our numerical calculations, unless otherwise indicated. This division of the transverse 
momentum phase space gives us practically the same results as the invariant mass phase 
space slicing technique. This was precisely checked by the lepton charge asymmetry results 
predicted by DYRAD p4|, and the NLO [up to 0{asy\ calculation within the ResBos Monte 
Carlo package. 

To facilitate our comparison, we calculate the NLO and the resummed distributions using 
the same parton luminosities and parton distribution functions, EW and QCD parameters, 
and renormalization and factorization scales so that any difference found in the distributions 
is clearly due to the different QCD physics included in the theoretical calculations. (Recall 
that they are different models of calculations based upon the same QCD theory, and the 
resummed calculation contains the dynamics of the multiple soft gluon radiation.) This way 
we compare the resummed and the NLO results on completely equal footing. The parton 
distributions used in the different order calculations are listed in Table In Table || and 
the rest of this work, we denote by Resummed 0{a'g) the result of the resummed calculation 
with and C(O'i) included [cf. Appendix D]; by Resummed 0{as) with A^^\ B^^^ 

and C(°'i); by Resummed O(a^) with A'^^\ B^^^ and (Similarly, later in Table |IT|, CSS 
0{ag) implies that A^^''^\ B^^''^^ and C^*^'^^ included in the resummation calculation, etc.) In 
the following, we discuss the relevant experimental observables predicted by these models of 
calculations using the ResBos code. Our numerical results are given for the Tevatron, a pp 
collider with \fS = 1.8 TeV, and CTEQ4 PDF's unless specified otherwise. 
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Fixed order 






Resummed 






0{oPs) 


0{as) 




0{oPs) 


0[as) 




PDF 


CTEQ4L 


CTEQ4M 


CTEQ4M 


CTEQ4L 


CTEQ4M 


CTEQ4M 



TABLE II. List of PDF's used at the different models of calculations. The values of the 
strong coupling constants used with the CTEQ4L and CTEQ4M PDF's are a^^^M^o) = 0.132 
and (M^o) = 0.116 respectively. 



A. Vector Boson Transverse Momentum Distribution 

According to the parton model the primordial transverse momenta of partons entering 
into the hard scattering are zero. This implies that a 7*, or boson produced in 
the Born level process has no transverse momentum, so that the LO Qt distribution is a 
Dirac-delta function peaking at Qt = 0. In order to have a vector boson produced with a 
non-zero Qt-, an additional parton has to be emitted from the initial state partons. This 
happens in the QCD process. However the singularity at Qt = prevails up to any fixed 
order in as of the perturbation theory, and the transverse momentum distribution da/dQ"^ 
is proportional to Q^^ x [1 or ln(Qf^/Q^)] at small enough transverse momenta. The most 
important feature of the transverse momentum resummation formalism is to correct this 
unphysical behavior and render da/dQ'^ finite at zero Qt- The Qt distributions of the 
and bosons predicted by the NLO and the resummed calculations are shown in Fig. 0. 

We find that in the resummed calculation, after taking out the resonance weighting factor 
Q'^ / {{Q"^ — My)^ + Q'^Ty / My) lu Eq. (|1|), the shape of the transverse momentum distribution 
of the vector boson V between Qt = and 20 GeV is remarkably constant for Q being in 
the vicinity of My. Fixing the rapidity y of the vector boson V at some value yo and taking 
the ratio 



R{Qt,Qo) = 



da 




dQ^dQldy 


Q=Qo,y=yo 


da 




dQ^dQldy 


Q=Mv,y=yo 



we obtain almost constant curves (within 3 percent) for Q = My ± 10 GeV (cf. Fig. ^ for 
V = and Z^. The fact that the shape of the transverse momentum distribution shows 
such a weak dependence on the invariant mass Q in the vicinity of the vector boson mass 
can be used to make the Monte Carlo implementation of the resummation calculation faster. 
This weak dependence was also used in the D0 W mass analysis when assuming that the 
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5 10 15 20 25 30 30 40 50 60 70 80 

Qj(GeV) Qj(GeV) 



FIG. 2. The low and intermediate Qt regions of the and distributions at the Tevatron, 

I) 



calculated in fixed order 0{as) (dotted) and 0(a|) (dash-dotted), and resummed 0{as) (dashed) 



and 0{a'g) (solid) [cf. Table III]. The cross-over occurs at 54 GeV for the 0{as)-, and at 49 GeV 
for the 0{al) distribut ions. The situation is very similar for the boson. 



mass dependence of the fully differential W boson production cross section factorizes as a 
multiplicative term Similarly, we define the ratio 



da 



R{QT,yo) 



dQ^dQldy 



=Mv,y=yo 



da 




dQ^dQldy 


Q=Mv,y=(} 



to study the Qt shape variation as a function of the vector boson rapidity. Our results are 
shown in Fig. ^. Unlike the ratio R{Qt,Qo) shown in Fig. |^, the distributions of R{QT,yo) 
for the and bosons are clearly different for any value of the rapidity yo. 

To utilize the information on the transverse momentum of the boson in Monte Carlo 
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FIG. 3. The ratio R{Qt,Qo), with uq = 0, for and bosons as a function of Qt- For 
W~^, sohd hnes are: Qq = 78 GeV (upper) and 82 GeV (lower), dashed: Qq = 76 GeV (upper) and 
84 GeV (lower), dotted: Qq = 70 GeV (upper) and 90 GeV (lower). For bosons, solid lines: 
Qo = 88 GeV (upper) and 92 GeV (lower), dashed: Qo = 86 GeV (upper) and 94 GeV (lower), 
dotted: Qo = 80 GeV (upper) and 100 GeV (lower). 



simulations to reconstruct the mass of the boson, it was suggested in Ref. to predict 
QriW^) distribution from the measured Qt{Z^) distribution and the calculated ratio of 
QriW^) and Qt{Z^) predicted by the resummation calculations [0,|12|5 in which the vector 
boson is assumed to be on its mass-shell. Unfortunately, this idea will not work with a good 
precision because, as clearly shown in Fig. ^, the ratio of the and transverse momen- 
tum distributions depends on the rapidities of the vector bosons. Since the rapidity of the 
boson cannot be accurately reconstructed without knowing the longitudinal momentum 
(along the beam pipe direction) of the neutrino, which is in the form of missing energy car- 
ried away by the neutrino, this dependence cannot be incorporated in data analysis and the 
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FIG. 4. The ratio R{QT,yo), with Qq = My, for and bosons as a function of Qt- 

above ansatz cannot be realized in practice for a precision measurement of Mi^/.^ Only the 
Monte Carlo implementation of the exact matrix element calculation (ResBos) can correctly 
predict the distributions of the decay leptons, such as the transverse mass of the W"^ boson, 
and the transverse momentum of the charged lepton, so that they can be directly compared 
with experimental data to extract the value of Mw- We comment on these results later in 
this section. 

Another way to compare the results of the resummed and the NLO calculations is given 
by the distributions of <7{Qt > Qt^^) / crTotah as shown in Fig. ^. We defined the ratio as 



^If a high precision measurement were not required, then one could choose from the two-fold 
solutions for the neutrino longitudinal momentum to calculate the longitudinal momentum of the 
boson. 
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^ _ a{QT > Qf^) 1 [QT"- ,^ da{h^h2 ^ V) 
Kcss = = / "Vr 



CTotal O-Total J QT' dQx 

where Q^^^ is the largest Qt allowed by the phase space. In the NLO calculation, a^Qx > 
Qmin-j gj-Q-^g without bound near Q™™ = 0, as the result of the singular behavior in the 

matrix element. The NLO curve runs well under the resummed one in the 2 GeV < Q™™ < 
30 GeV region, and the Qt distributions from the NLO and the resummed calculations have 
different shapes even in the region where Qt is of the order 15 GeV. 

With large number of fully reconstructed events at the Tevatron, one should be able 
to use data to discriminate these two theory calculations. In view of this result it is not 
surprising that the D0 analysis of the as measurement based on the measurement 
of a{W + 1 jet)/a{W + jet) does not support the NLO calculation in which the effects 
of the multiple gluon radiation are not included. We expect that if this measurement was 
performed by demanding the transverse momentum of the jet to be larger than about 50 
GeV, at which scale the resummed and the NLO distributions in Fig. |] cross, the NLO 
calculation would adequately describe the data. 

To show that for Qt below 30 GeV, the QCD multiple soft gluon radiation is important 
to explain the D0 data p6|, we also include in Fig. |^ the prediction for the Qt distribution 
at the order of a|. As shown in the figure, the q;| curve is closer to the resummed curve 
which proves that for this range of Qt the soft gluon effect included in the a| calculation is 
important for predicting the vector boson Qt distribution. In other words, in this range of 
Qt, it is more likely that soft gluons accompany the boson than just a single hard jet 
associated with the vector boson production. For large Qt, it becomes more likely to have 
hard jet(s) produced with the vector boson. 

Measuring Rcss in the low Qt region (for Qt ^ Q/2) provides a stringent test of the 
dynamics of the multiple soft gluon radiation predicted by the QCD theory. The same 
measurement of Rcss can also provide information about some part of the non-perturbative 
physics associated with the initial state hadrons. As shown in Fig. |] and in Ref. p2|], the 
effect of the non-perturbative physics on the Qt distributions of the and Z° bosons 
produced at the Tevatron is important for Qt less than about 10 GeV. This is evident by 
observing that different parametrizations of the non-perturbative functions do not change 
the Qt distribution for Qt > 10 GeV, although they do dramatically change the shape of 
Qt for Qt < 10 GeV. Since for and Z° production, the \ia{Q'^/Ql) term is large, the 
non-perturbative function, as defined in Eq. (^, is dominated by the Fi{b) term (or the g2 
parameter) which is supposed to be universal for all Drell-Yan type processes and related to 
the physics of the renormalon ||20|. Hence, the measurement of Rcss cannot only be used to 
test the dynamics of the QCD multiple soft gluon radiation, in the 10 GeV < Qt < 40 GeV 
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FIG. 5. The ratio Rcss as a function of Q™" for and bosons. The fixed order {0{as) 
short dashed, 0{a'g) dashed) curves are ill-defined in the low Qt region. The resummed (solid) 
curves are calculated for 52 = 0.38 (low), 0.58 (middle) and 0.68 (high) GeV^ values. 

region, but may also be used to probe this part of non-perturbative physics for Qt < 10 
GeV. It is therefore important to measure Rcss at the Tevatron. With a large sample of 
data at Run 2, it is possible to determine the dominant non-perturbative function which 
can then be used to calculate the boson Qt distribution to improve the accuracy of the 
Mw and the charged lepton rapidity asymmetry measurements. 

B. The Total Cross Section 



Before we compare the distributions of the decay leptons, we examine the question 
whether or not the Qt resummation formalism changes the prediction for the total cross 
section. In Ref. it was shown that in the AEGM formalism, which differs from the CSS 
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FIG. 6. Transverse momentum distributions of and bosons calculated with low (long 
dash, 52 = 0.38 GeV^), nominal (solid, 52 = 0.58 GeV^) and high (short dash, 52 = 0.68 GeV^) 52 
non-perturbative parameter values. The low and high excursions in are the present one standard 
deviations from the nominal value in the Ladinsky-Yuan parametrization. 

formalism, the 0{as) total cross section is obtained after integrating their resummation 
formula over the whole range of the phase space. 

In the CSS formalism, without including the C and Y functions, the fully integrated 
resummed result recovers the 0{a%) cross section, provided that Qt is integrated from zero 
to Q. This can be easily verified by expanding the resummation formula up to 0{as).i 
dropping the C^^^ and the Y pieces (which are of order 0{as)), and integrating over the 
lepton variables. It yields 




TT 



B^'^ In ( S 
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(^siQ) (Q 



27r -"\P? 



T/ 

+ fj/hA^i. Q') {P-k^t ® h/h,) ix2, Q')] +3^k}, (9) 

where Pt is the upper hmit of the Qt integral and we fixed the mass of the vector boson for 
simphcity. To derive the above resuh we have used the canonical set of the Ci {i = 1,2,3) 
coefficients (cf. Appendix C and D). When the upper limit Pt is taken to be Q, all the logs 
in the above equation vanish and Eq. reproduces the Born level {0{a^)) cross section. 
Similar conclusion holds for higher order terms from the expansion of the resummed piece 
when C and Y are not included. This is evident because the singular pieces from the 
expansion are given by 



da 



1 oo 2n-l ' 
VT n=^ m=n \ V , 



singular n=l m=0 

The integral of these singular terms will be proportional to \n{Q'^/P^) raised to some power. 
Again, for Pt = Q all the logs vanish and the tree level result is obtained. 

The inclusion of C-^-* and Y functions will change the above conclusion and lead to a 
different total cross section, because C*^^-' contains the hard part virtual corrections and Y 
contains the hard gluon radiation. In the Qt < Q region the resummed piece dominates, 
because it resums the singular pieces: x [1 or ln(Q^/(5|,)], to all order. However, in 
the Qt > Q region the perturbative result should be used, because the singular pieces do 
not dominate (the logs are small in this case) and the other contributions in the fixed order 
perturbative contributions can be important. In order to define the total cross section a 
prescription of smooth transition between the resummed and the fixed order perturbative 
results is necessary. In Fig. H, we show the resummed (1,1) (resummation with A^^\ B^^^ and 
Qio,i) included) and the fixed order 0{as) Qt distributions for and bosons. As shown, 
the resummed (1,1) and the fixed order curves are close to each other for Q/2 < Qt < Q, 
and they cross near Qt ~ Q/2. Henceforth, we adopt the following matching procedure. 
For Qt values below the crossing point we use the resummed curve, and above it the 0{as) 
curve, to define a resummed 0{as) distribution. This distribution is smooth (although 
not differentiable at the matching point), and most importantly it does not alter either the 
resummed or the 0{as) distributions in the kinematic regions where they are proven to be 
valid. This simple matching prescription provides us with an 0{as) resummed total cross 
section with an error of 0{a1), as shown in Ref. |]12|. In practice this translates into less than 
a percent deviation between the fixed order and the resummed 0{as) total cross sections. 
This can be understood from the earlier discussion that if the matching were done at Qt 
equal to Q, then the total cross section calculated from the CSS resummation formalism 
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11.3 
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w+ 


2.0 


9.71 


12.5 


12.6 


12.6 


12.7 






1.8 


5.23 


6.69 


6.79 


6.79 


6.82 


6.86 ± 0.36 




2.0 


6.11 


7.47 


7.52 


7.52 


7.57 





TABLE III. Total cross sections of pp — > (W~^or Z^)X at the present and upgraded Tevatron, 
calculated in different prescriptions, in units of nb. The finite order total cross section results are 
based on the calculations in Ref. [27]. The Pert. 0(a|) results were obtained from Ref. [28]. The 
"+" signs refer to the inclusion of the Y piece and the "©" signs to the switch from the resummed 
to the fixed order calculations. 



should be the same as that predicted by the NLO calculation, provided that C*-^-* and Y^'^^ 
are included. However, this matching prescription would not result in a smooth curve for 
the Qt distribution at Qt = Q- The small difference between the resummed 0{as) and the 
NLO [fixed order, 0{as)] total cross sections comes from the matching procedure described 
above. This difference indicates the size of the higher order corrections not included in the 
NLO calculation. 

In order to compare the Qt distribution with experimental data we also include the effect 
of some known higher order corrections to the Sudakov factor and plot the resummed (2,1) 
(with and C(O'i) included) and the fixed order 0(a|) distributions [||].0 We 

match the two distributions at the value of Qt, about 48 GeV, where they cross over for 
production, that defines the resummed 0(a|) distribution in the whole Qt region. The 
total cross section predicted from various theory calculations are listed in Table 

We note that kinematic cuts affect the total cross section in a subtle manner. It is obvious 
from our matching prescription that the resummed 0(a|) and the fixed order 0{as) curves 
in Fig. 1^ will never cross. On the other hand, the resummed 0(a|) total cross section is about 
the same as the fixed order 0{as) cross section when integrating Qt from to Q. These 
two facts imply that when kinematic cuts are made on the Qt distribution with Qt < Q, 
we will obtain a higher total cross section in the fixed order 0{as) than in the resummed 



^The K factor, which is defined to be the ratio of the cross sections at 0(a|) and 0{as), as the 
function of Qt, changes within about 5% in the plotted Qt region. For W~^: K{Qt = 30 GeV) = 
1.48, K(50 GeV) = 1.43, K(70 GeV) = 1.39; and for Z^: K{Qt = 30 GeV) = 1.47, K{50 GeV) = 
1.43, K(70 GeV) = 1.39. 
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C(a;|) calculation. In this paper we follow the CDF cuts (for the boson mass analysis) 
and demand Qt < 30 GeV Consequently, in many of our figures, to be shown below. 



the fixed order 0{as) curves give about 3% higher total cross section than the resummed 
ones. 

C. Lepton Charge Asymmetry 



The CDF lepton charge asymmetry measurement ||29| played a crucial role in constraining 



the slope of the u/d ratio in recent parton distribution functions. It was shown that one 
of the largest theoretical uncertainty in the mass measurement comes from the parton 
distributions [^, and the lepton charge asymmetry was shown to be correlated with the 
transverse mass distribution [^]. Among others, the lepton charge asymmetry is studied to 
decrease the errors in the measurement of Mw coming from the parton distributions. Here 
we investigate the effect of the resummation on the lepton rapidity distribution, although it 
is not one of those observables which are most sensitive to the Qt resummation, i.e. to the 
effect of multiple soft gluon radiation. 

The definition of the charge asymmetry is 

da da 



+ 



dy+ dy_ 

where y+ (?/_) is the rapidity of the positively (negatively) charged particle (either vector 
boson or decay lepton). Assuming CP invariance,[] the following relation holds: 

da da 

[y) = -j—{-y)- 



dy+ dy_ 
Hence, the charge asymmetry is frequently written as 

da , . da , 

dcr, ^, da , ^, 
-r^y > 0) + -r{y < 0) 
dy dy 

For the charge asymmetry of the vector boson (VT^) or the charged decay lepton (i^), the 
fixed order and the resummed 0{as) [or (9(a|)] results are the same, provided that there 



Here we ignore the small CP violating effect due to the CKM matrix elements in the SM. 
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FIG. 7. Lepton charge asymmetry distributions, (a) Without any kinematic cuts, the NLO 
(long dashed) and the resummed 0{as) (sohd) curves overlap and the LO (short dashed) curve 
differs somewhat from them, (b) With cuts {Qt < 30 GeV,p^ ''^ > 25 GeV), the effect of the 
different Qt distributions renders the lepton rapidity asymmetry distributions different. The two 
resummed curves calculated with §2 = 0.58 and 0.78 GeV^^ cannot be distinguished on this plot. 

are no kinematic cuts imposed. This is because the shape difference in the vector boson 
transverse momentum has been integrated out and the total cross sections are the same 
up to higher order corrections in (as)- In Fig. |^(a) we show the lepton charge asymmetry 
without cuts for CTEQ4M PDF. The NLO and the resummed curves overlap, although they 
differ from the 0{ag) prediction. 

On the other hand, when kinematic cuts are applied to the decay leptons, the rapidity 
distributions of the vector bosons or the leptons in the fixed order and the resummed cal- 
culations are different. Restriction of the phase space implies that only part of the vector 
boson transverse momentum distribution is sampled. The difference in the resummed and 
the fixed order Qt distributions will prevail as a difference in the rapidity distributions of 
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the charged lepton. We can view this phenomenon in a different (a Monte Carlo) way. In the 
rest frame of the W^, the decay kinematics is the same, whether it is calculated up to 0{as) 
or within the resummation formalism. On the other hand, the rest frame is different for 
each individual Monte Carlo (MC) event depending on the order of the calculation. This dif- 
ference is caused by the fact that the Qt distribution of the is different in the 0{as) and 
the resummed calculations, and the kinematic cuts select different events in these two cases. 
Hence, even though the Qt distribution of the is integrated out, when calculating the 
lepton rapidity distribution, we obtain slightly different predictions in the two calculations. 
The difference is larger for larger \y^\, being closer to the edge of the phase space, because 
the soft gluon radiation gives high corrections there and this effect up to all order in as is 
contained in the resummed but only up to order of as in the NLO calculation. Because the 
rapidity of the lepton and that of the vector boson are highly correlated, large rapidity lep- 
tons mostly come from large rapidity vector bosons. Also, a vector boson with large rapidity 
tends to have low transverse momentum, because the available phase space is limited to low 
Qt for a boson with large \y\. Hence, the difference in the low Qt distributions of the 
NLO and the resummed calculations yields the difference in the distribution for leptons 
with high rapidities. 

Asymmetry distributions of the charged lepton with cuts using the CTEQ4M PDF are 
shown in Fig. 0(b). The applied kinematic cuts are: Qt < 30 GeV, Pt > 25. These are 
the cuts that CDF used when extracted the lepton rapidity distribution from their data 



We have checked that the ResBos fixed order 0{as) curve agrees well with the DYRAD 
result. As anticipated, the 0{ag), 0{as) and resummed results deviate at higher rapidities 
{{y'^l > 1.5).0 The deviation between the NLO and the resummed curves indicates that to 
extract information on the PDF in the large rapidity region, the resummed calculation, in 
principle, has to be used if the precision of the data is high enough to distinguish these 
predictions. Fig. ^(b) also shows the neghgible dependence of the resummed curves on the 
non-perturbative parameter g2. We plot the result of the resummed calculations with the 
nominal g2 = 0.58 GeV~^, and with g2 = 0.78 GeV~^ which is two standard deviations 
higher. The deviation between these two curves (which is hardly observable on the figure) 
is much smaller than the deviation between the resummed and the NLO ones. 

There is yet another reason why the lepton charge asymmetry can be reliably predicted 
only by the resummed calculation. When calculating the lepton distributions in a numerical 



°Here and henceforth, unless specified otherwise, by a resummed calculation we mean our re- 
summed 0{a'g) result. 
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FIG. 8. Distributions of positron rapidities from the decays of VF"'"'s produced at the Tevatron, 
predicted by the resummed (sohd) and the NLO (dashed) calculations with the same kinematic 
cuts as for the asymmetry plot. 

0{as) code, one has to artificially divide the vector boson phase space into hard and soft 
regions, depending on - for example - the energy or the Qt of the emitted gluon (e.g. 
qq — > W+ hard or soft gluon). The observables calculated with this phase space slicing 
technique acquire a dependence on the scale which separates the hard from the soft regions. 
To emphasize this dependence as an example, we show that when the phase space is divided 
by the Qt separation, the dependence of the asymmetry on the scale Q^^^ can be comparable 
to the difference in the 0{as) and the resummed results. This means that there is no definite 
prediction from the NLO calculation for the lepton rapidity distribution. Only the resummed 
calculation can give an unambiguous prediction for the lepton charge asymmetry. 

Before closing this section, we also note that although in the lepton asymmetry distribu- 
tion the NLO and resummed results are about the same for \y'^'^\ < 1, it does not imply that 
the rapidity distributions of the leptons predicted by those two theory models are the same. 
As shown in Fig. |, this difference can in principle be observable with a large statistics data 
sample and a good knowledge of the luminosity of the colliding beams. 



Since the invariant mass of the boson cannot be reconstructed without knowing the 
longitudinal momentum of the neutrino, one has to find a quantity that allows an indirect 
determination of the mass of the boson. In the discovery stage of the bosons at the 
CERN SppS collider, the mass and width were measured using the transverse mass distri- 



D. Transverse Mass Distribution 
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bution of the charged lepton-neutrino pair from the boson decay. Ever since the early 



eighties, the transverse mass distribution, m-p = ^J2pj,pj,{l — cos A0ei/), has been known as 
the best measurable for the extraction of both Mw and Tw, for it is insensitive to the trans- 
verse momentum of the boson. The effect of the non-vanishing vector boson transverse 



momentum on the distribution was analyzed p2| , p3| well before the Qt distribution 
of the boson was correctly calculated by taking into account the multiple soft gluon 



radiation. Giving an average transverse boost to the vector boson, the authors of Ref. |^| 
concluded that for the fictive case of T\y = 0, the end points of the transverse mass distri- 
bution are fixed at zero: da / drn^{m'^ = 0) = da/dm^{'m\. = M^) = 0. The sensitivity of 
the rriT shape to a non-zero Qt is in the order of {{Qt/MwY) ~1% without affecting the 
end points of the rriT distribution. Including the effect of the finite width of the boson, 
the authors in Ref. ||33| showed that the shape and the location of the Jacobian peak are 
not sensitive to the Qt of the boson either. The non-vanishing transverse momentum 
of the boson only significantly modifies the rriT distribution around rriT = 0. 

Our results confirm that the shape of the Jacobian peak is quite insensitive to the order of 
the calculation. We show the NLO and the resummed transverse mass distributions in Fig. |^ 
for bosons produced at the Tevatron with the kinematic cuts: Qt < 30 GeV, > 25 
GeV and \y^^\ < 3.0. Fig. ^(a) covers the full (experimentally interesting) ttit range while 
Fig. ^(b) focuses on the ttlt range which contains most of the information about the 
mass. There is little visible difference between the shapes of the NLO and the resummed 
rriT distributions. On the other hand, the right shoulder of the curve appears to be "shifted" 
by about 50 MeV, because, as noted in Section [111 B[ the total cross sections are different 
after the above cuts imposed in the NLO and the resummed calculations. At Run 2 of the 
Tevatron, with large integrated luminosity (~ 2 fb~^), the goal is to extract the boson 



mass with a precision of 30-50 MeV from the rriT distribution |^ . Since Mw is sensitive to 
the position of the Jacobian peak [^, the high precision measurement of the mass has 
to rely on the resummed calculations. 

Extraction of Mw from the transverse mass distribution has some drawbacks though. The 
reconstruction of the transverse momentum p^ of the neutrino involves the measurement of 
the underlying event transverse momenta: p^^ = — j>/ — — jj^^^.'^^'-v^^a event ^ This 

resolution degrades by the number of interactions per crossing (A''/J With a high 

luminosity (~ 100 fb"^) at the 2 TeV Tevatron (TEV33) |4[, Nj^ can be as large as 10, 
so that the Jacobian peak is badly smeared. This will lead to a large uncertainty in the 
measurement of Mw For this reason the systematic precision of the reconstruction will 
be less at the high luminosity Tevatron, and an Mw measurement that relies on the lepton 
transverse momentum distribution alone could be more promising. We discuss this further 
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FIG. 9. Transverse mass distribution for production and decay at the 1.8 TeV Tevatron. 



in the next section. 

The theoretical Hmitation on the measurement using the ttit distribution comes from 
the dependence on the non-perturbative sector, i.e. from the PDF's and the non-perturbative 
parameters in the resummed formahsm. Assuming the PDF's and these non-perturbative 
parameters to be independent variables, the uncertainties introduced are estimated to be 



less than 50 MeV and 10 MeV, respectively, at the TEV33 1^ is clear that the 

main theoretical uncertainty comes from the PDF's. As to the uncertainty due to the non- 
perturbative parameters (e.g. (72) in the CSS resummation formalism, it can be greatly 
reduced by carefully study the Qt distribution of the boson which is expected to be 
copiously produced at Run 2 and beyond. 

The Mvi/ measurement at the LHC may also be promising. Both ATLAS and CMS 
detectors are well optimized for measuring the leptons and the missing The cross 

section of the boson production is about four times larger than that at the Tevatron, 
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and in one year of running with 20 fb~^ luminosity yields a few times 10^ W ^ iu events 
after imposing similar cuts to those made at the Tevatron. Since the number of interactions 
per crossing may be significantly lower (in average Nj^ = 2) at the same or higher luminosity 
than that at the TEV33 [0, the Jacobian peak in the rriT distribution will be less smeared 
at the LHC than at the TEV33. Furthermore, the non-perturbative effects are relatively 
smaller at the LHC because the perturbative Sudakov factor dominates. On the other hand, 
the probed region of the PDF's at the LHC has a lower value of the average x (~ 10~^) than 
that at the Tevatron (~ 10'^), hence the uncertainty from the PDF's might be somewhat 
larger. A more detailed study of this subject is desirable. 

E. Lepton Transverse Momentum 

Due to the limitations mentioned above, the transverse mass (my) method may not be 
the only and the most promising way for the precision measurement of Mw at some future 
hadron colliders. As discussed above, the observable niT was used because of its insensitivity 
to the high order QCD corrections. In contrast, the lepton transverse momentum (p^) 
distribution receives a large, O {{Qt / My/)) ~ 10%perturbative QCD correction at the order 
0:5, as compared to the Born process. With the resummed results in hand it becomes possible 
to calculate the pf. distribution precisely within the perturbative framework, and to extract 
the mass straightly from the transverse momentum distributions of the decay leptons. 

Just like in the itit distribution, the mass of the boson is mainly determined by the 
shape of the distribution near the Jacobian peak. The location of the maximum of the peak 
is directly related to the boson mass, while the theoretical width of the peak varies with 
its decay width V-^. Since the Jacobian peak is modified by effects of both Qt and it 
is important to take into account both of these effects correctly. In our calculation (and in 
ResBos) we have properly included both effects. 

The effect of resummation on the transverse momentum distribution of the charged lepton 
from and decays is shown in Fig. [TO. The NLO and the resummed distributions differ 



a great amount even without imposing any kinematic cuts. The clear and sharp Jacobian 
peak of the NLO distribution is strongly smeared by the finite transverse momentum of 
the vector boson introduced by multiple gluon radiation. This higher order effect cannot 
be correctly calculated in any finite order of the perturbation theory and the resummation 
formalism has to be used. 

One of the advantages of using the distribution to determine My/ is that there is 
no need to reconstruct the distribution which potentially limits the precision of the ttit 
method. From the theoretical side, the limitation is in the knowledge of the non-perturbative 
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FIG. 10. Transverse momentum distributions of from and decays for the NLO 
(dashed) and the resummed 0{as) (sohd) calculations. Resumming the initial state multiple 
soft-gluon emission has the typical effect of smoothening and broadening the Jacobian peak (at 

= My/2). The CDF cuts are imposed on the distributions, but there are no cuts on the 

distributions. 

sector. Studies at D0 show that the distribution is most sensitive to the PDF's 
and the value of the non-perturbative parameter g2. The distribution is more sensitive 
to the PDF choice, than the rriT distribution is. The uncertainty in the PDF causes an 
uncertainty in M]^ of about 150 MeV, which is about three times as large as that using the 
rriT method A 0.1 GeV^ uncertainty in g2 leads to about AM^ = 30 MeV uncertainty 
from the p^ fit, which is about five times worse than that from the itlt measurement ||35|| . 
Therefore, to improve the Mw measurement, it is necessary to include the data sample at 
the high luminosity Tevatron to refit the Qi^s and obtain a tighter constrain on them from the 
Qt distribution of the Z^ boson. The D0 study showed that an accuracy of Ag2 = 0.01 GeV^ 
can be achieved with Run 2 and TeV33 data, which would contribute an error of AMw < 5 
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MeV from the pf, |3^. In this case the uncertainty coming from the PDF's remains to be the 
major theoretical hmitation. At the LHC, the distribution can be predicted with an even 
smaller theoretical error coming from the non-perturbative part, because at higher energies 
the perturbative Sudakov factor dominates over the non-perturbative function. 

It was recently suggested to extract from the ratios of the transverse momenta of 
leptons produced in and Z° decay [^. The theoretical advantage is that the non- 
perturbative uncertainties are decreased in such a ratio. On the other hand, it is not enough 
that the ratio of cross sections is calculated with small theoretical errors. For a preci- 
sion extraction of the mass the theoretical calculation must be capable of reproducing 
the individually observed transverse momentum distributions themselves. The mass 
measurement requires a detailed event modeling, understanding of detector resolution, kine- 
matical acceptance and efficiency effects, which are different for the and events, as 
illustrated above. Therefore, the ratio of cross sections can only provide a useful check for 
the mass measurement. 

For Drell-Yan events or lepton pairs from Z'^ decays, additional measurable quantities 
can be constructed from the lepton transverse momenta. They are the distributions in the 
balance of the transverse momenta {Apr = \Pt\ ~ \Pt\) ^^e angular correlation of the 
two lepton momenta {z = —p^ ■ p^ /\max.{p^ ,p^y\^). It is expected that these quantities 
are also sensitive to the effects of the multiple soft gluon radiation. These distributions are 



shown in Figure |Tl|. As shown, the resummed distributions significantly differ from the NLO 
ones. In these, and the following figures for Z° decay distributions, it is understood that the 
following kinematic cuts are imposed: Qf" < 30 GeV, > 25 GeV and | < 3.0, 

unless indicated otherwise. 



F. Lepton Angular Correlations 

Another observable that can serve to test the QCD theory beyond the fixed-order per- 
turbative calculation is the difference in the azimuthal angles of the leptons and [2 from 
the decay of a vector boson V . In practice, this can be measured for 7* or Z^ ^1^2- We 
show in Fig. the difference in the azimuthal angles of e"*" and e~ [A.(j)^^^ ), measured in 
the laboratory frame for Z^ — > e"'"e~, calculated in the NLO and the resummed approaches. 
As indicated, the NLO result is ill-defined in the vicinity of A0 ~ vr, where the multiple 
soft-gluon radiation has to be resummed to obtain physical predictions. 

Another interesting angular variable is the lepton polar angle distribution cos^*^ in the 
Collins-Soper frame. It can be calculated for the Z^ decay and used to extract sin at 
the Tevatron The asymmetry in the polar angle distribution is essentially the same 



29 




FIG. 11. Balance in transverse momentum Apj- = \prp-\ — \Px\ and angular correlation 
z = —p^ • p^^/[max(p^ of the decay leptons from bosons produced at the Tevatron. 



as the forward-backward asymmetry Aps measured at LEP. Since Aps depends on the 
invariant mass Q and around the energy of the peak ApB happens to be very small, the 
measurement is quite challenging. At the hadron collider, on the other hand, the invariant 



mass of the incoming partons is distributed over a range so the asymmetry is enhanced [15 



The potentials of the measurement deserve a separated study. In Fig. ^ we show the 
distributions of cos^*^ predicted from the NLO and the resummed results. 



G. Vector Boson Longitudinal Distributions 



The resummation of the logs involving the transverse momentum of the vector boson does 
not directly affect the shape of the longitudinal distributions of the vector bosons. A good 
example of this is the distribution of the longitudinal momentum of the Z*^ boson which can 
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FIG. 12. The correlation between the lepton azimuthal angles near the region Ac/) ~ vr for 
pp — > {Z^ — > e'^e~)X. The resummed (solid) distribution gives the correct angular correlation 
of the lepton pair. The NLO (dashed lines) distribution near Ac/) = tt is ill-defined and depends 
on (the scale for separating soft and hard gluons in the NLO calculation). The two NLO 

distributions were calculated with Q^^^^ = 1.2 GeV (long dash) and Q^'^ = 2.0 GeV (short dash). 

be measured at the Tevatron with high precision, and can be used to extract information 
on the parton distributions. It is customary to plot the rescaled quantity Xp = 2q^/\/S, 
where = smh{y) Q'^ + Qj, is the longitudinal momentum of the boson measured in 
the laboratory frame. In Fig. |14|, we plot the distributions predicted in the resummed and 
the NLO calculations. As shown, their total event rates are different in the presence of 
kinematic cuts. (Although they are the same if no kinematic cuts imposed.) This conclusion 
is similar to that of the distributions, as discussed in Sections [III B| and [III C . 

Without any kinematic cuts, the vector boson rapidity distributions are also the same 
in the resummed and the NLO calculations. This is so because when calculating the y 
distribution the transverse momentum Qt is integrated out so that the integral has the 
same value in the NLO and the resummed calculations. On the other hand, experimental 
cuts on the final state leptons restrict the phase space, so the difference between the NLO 
and the resummed Qt distributions affects the vector boson rapidity distributions. This 
shape difference is very small at the vector boson level, as shown in Fig. |15[ 

IV. DISCUSSION AND CONCLUSIONS 

With a 100 pb~^ luminosity at the Tevatron, around 2 x 10^ and 6 x 10^ bosons 
are produced, and the data sample will increase by a factor of 20 in the Run 2 era. In view 
of this large event rate, a careful study of the distributions of leptons from the decay of the 
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FIG. 13. Distribution of the e"*" polar angle cos{9*) in the Collins-Soper frame from decays 
at the Tevatron with cuts indicated in the text. 

vector bosons can provide a stringent test of the rich dynamics of the multiple soft gluon 
emission predicted by the QCD theory. Since an accurate determination of the mass of the 
boson and the test of parton distribution functions demand a highly precise knowledge 
of the kinematical acceptance and the detection efficiency of or bosons, the effects of 
the multiple gluon radiation have to be taken into account. In this work, we have extended 
the formalism introduced by Collins, Soper and Sterman for calculating an on-shell vector 
boson to include the effects of the polarization and the decay width of the vector boson 
on the distributions of the decay leptons. Our resummation formalism can be applied to 
any vector boson V where = 7*, , Z^, W, Z', etc., with either vector or axial- vector 
couplings to fermions (leptons or quarks). To illustrate how the multiple gluon radiation 
can affect the distributions of the decay leptons, we studied in detail various distributions 
for the production and the decay of the vector bosons at the Tevatron. 

One of the methods to test the rich dynamics of the multiple soft gluon radiation predicted 
by the QCD theory is to measure the ratio Rcss = '^^^l^^^ ^ for the and Z^ bosons. We 
found that, for the vector boson transverse momentum less than about 30 GeV, the difference 
between the resummed and the fixed order predictions (either at the as or a| order) can be 
distinguished by experimental data. This suggests that in this kinematic region, the effects 
of the multiple soft gluon radiation are important, hence, the Qt distribution of the vector 
boson provides an ideal opportunity to test this aspect of the QCD dynamics. For Qt less 
than about 10 GeV, the distribution of Qt is largely determined by the non-perturbative 
sector of QCD. At the Tevatron this non-perturbative physics, when parametrized by Eq. (^ 
for and Z^ production, is dominated by the parameter (72 which was shown to be related 
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FIG. 14. Longitudinal xp distributions of bosons produced at the Tevatron. The NLO 
(dashed) curves overestimate the rate compared to the resummed (sohd) ones, because kinematic 
cuts enhance the low Qt region where the NLO and resummed distributions are qualitatively 
different. Without cuts, the NLO and the resummed xp distributions are the same. 



to properties of the QCD vacuum |^0[. Therefore, precisely measuring the Qt distribution 
of the vector boson in the low Qt region, e.g. from the ample events, can advance our 
knowledge of the non-perturbative QCD physics. 

Although the rapidity distributions of the leptons are not directly related to the transverse 
momentum of the vector boson, they are predicted to be different in the resummed and the 
fixed order calculations. This is because to compare the theoretical predictions with the 
experimental data, some kinematic cuts have to be imposed so that the signal events can 
be observed over the backgrounds. We showed that the difference is the largest when the 
rapidity of the lepton is near the boundary of the phase space (i.e. in the large rapidity 
region), and the difference diminishes when no kinematic cuts are imposed. When kinematic 
cuts are imposed another important difference between the results of the resummed and the 
NLO calculations is the prediction of the event rate. These two calculations predict different 
normalizations of various distributions. For example, the rapidity distributions of charged 
leptons {y^^) from the decays of bosons are different. They even differ in the central 
rapidity region in which the lepton charge asymmetry distributions are about the same (cf. 
Figs. 1^ and ||). As noted in Ref. |^, with kinematic cuts, the measurement of Mw is 
correlated to that of the rapidity and its asymmetry through the transverse momentum of 
the decay lepton. Since the resummed and the NLO results are different and the former 
includes the multiple soft gluon emission dynamics, the resummed calculation should be 
used for a precision measurement of My/- 
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FIG. 15. Rapidity distributions (resummed: solid, NLO: dashed) of bosons produced at the 
Tevatron with the kinematic cuts given in the text. 

In addition to the rapidity distribution, we have also shown various distributions of the 
leptons which are either directly or indirectly related to the transverse momentum of the 
vector boson. For those which are directly related to the transverse momentum of the vector 
boson, such as the transverse momentum of the lepton and the azimuthal correlation of 
the leptons, our resummation formalism predicts significant differences from the fixed order 
perturbation calculations in some kinematic regions. The details were discussed in Section 
III. 

As noted in the Introduction, a full event generator, such as ISAJET, can predict a 
reasonable shape for various distributions because it contains the backward radiation algo- 
rithm p3|, which effectively includes part of the Sudakov factor, i.e. effects of the multiple 



gluon radiation. However, the total event rate predicted by the full event generator is usually 
only accurate at the tree level, as the short distance part of the virtual corrections cannot yet 
be consistently implemented in this type of Monte Carlo program. To illustrate the effects 
of the high order corrections coming from the virtual corrections, which contribute to the 



Wilson coefficients C in our resummation formalism, we showed in Fig. ^ the predicted 
distributions of the transverse momentum of the Drell-Yan pairs by ISAJET and by Res- 
Bos (our resummed calculation). In this figure we have rescaled the ISAJET prediction to 
have the same total rate as the ResBos result, so that the shape of the distributions can be 
directly compared. We restrict the invariant mass of the virtual photons Q to be between 
30 and 60 GeV without any kinematic cuts on the leptons. If additional kinematic cuts 
on the leptons are imposed, then the difference is expected to be enhanced, as discussed in 
Section [III Q . As clearly shown, with a large data sample in the future, it will be possible to 
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FIG. 16. Transverse momentum distribution of virtual photons in pp ^ ^* ^ e^e events 
predicted by ResBos (solid curve) and IS A JET (histogram), calculated for the invariant mass range 
30 GeV < Q < 60 GeV at the 1.8 TeV Tevatron. 

experimentally distinguish between these two predictions, and, more interestingly, to start 
probing the non-perturbative sector of the QCD physics. 
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APPENDIX A: KINEMATICS 



Here we summarize some details of the kinematics for the lepton pair production pro- 
cess hih2 iii2)X. The laboratory (lab) frame is the center-of-mass frame of the 
colliding hadrons hi and In the lab frame, the cartesian coordinates of the hadrons 
are: Pf^_^ /ijl^ct^) = (1, 0, 0, ±1), where is the center-of-mass energy of the collider. 
Transverse momentum resummation is performed in the Collins-Soper [CS) frame This 
is the special rest frame of the vector boson in which the z axis bisects the angle between 
the hi hadron momentum ph^CS) and the negative /12 hadron momentum —ph^{CS) ||19|| . 

To derive the Lorentz transformation A'^^(/a6 ^ CS) that connects the lab and CS 
frames (in the active view point): p^{CS) = h.^^{lab CS) p'^{lab), we follow the definition 
of the CS frame. First, we find the rotation which makes the azimuthal angle 0y of the 
vector boson vanish. Then, we find the boost into a vector boson rest frame. Finally, in the 
vector boson rest frame we find the rotation which brings the hadron momentum ph^{CS) 
and negative hadron momentum —ph2{CS) into the desired directions. 

The transformation that takes g'^(/a6) = (g'^, cos </)y, sin 0y, g'^) from the lab frame 
into a longitudinally co-moving frame (long), in which q^{long) = Qt, 0, q^), is a rotation 
around the z axis, which is 



M'^{lab long) 



(\ 

COS^y sin0y 

-Sin0y COS0y 

VO 1/ 



A boost by /3 = —q{long)j(f brings four vectors from the longitudinal [long) frame into a 
vector boson rest frame (rest) . The matrix of the Lorentz boost from the long frame to the 
rest frame, expressed explicitly in terms of q^ is 



A^j^{long rest) 



( 



\ 



Q 

Or 

Q 



~Q 



Or 

Q 

Ml 



Qiq' + Q) 



Qtq' 



^ 


1 









Qiq' + Q) 



Qtq' 



Q 



\ 



Qiq' + Q) 



Q{q' + Q) + {q - 
Q{q' + Q) 





3\2 



where Q = ^ (g°)^ — Qt ~ (q^)"^ is the vector boson invariant mass, and the transverse 
mass is defined as Mt = ^JQ"^ + Qt- The transformation from the lab frame to the rest 
frame is then the product of the above boost and rotation: K^^{lab —>■ rest) = A'^^{long —>■ 
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rest)A\{lab — > long). Had we used only the boost /3' = —q{lab)/q^, we would have obtained 



the same 



result for k^^ilah — > rest) in one step: 



k^^{lah rest) = 



-— sm0^ 

■2 



—sin 0y cos 0y 

,3 



g3 ^^ 



■ sm ( 



Q{q' + Q) 




After boosting the lab frame hadron momenta into this rest frame, we obtain 



Q{q' + Q) I 



Ph.Mi^^'t) = Kil^^b ^ rest) p^^^^ilab) = 



Q 



Q±q' (±g-a(g° + g)±(g-y^ 



g(g° + Q) 



and the polar angles of p^^(rest) and —p^^{rest) are not equal unless Qt = 0. (In the above 
expressions the upper signs refers to hi and the lower signs to /i2-) In the general Qt 7^ case 
we have to apply an additional rotation in the rest frame so that the z-axis bisects the angle 
between the hadron momentum phi{CS) and the negative hadron momentum ~ph2{CS). It 
is easy to see that to keep PhiM in the xz plane, this rotation should be a rotation around 
the y axis by an angle a = arccos [Q{q^ + M^) / {Mxiq^ + Q))]. 

Thus the Lorentz transformation from the lab frame to the CS frame is A^j,(/a6 ^ CS) = 
h^y^{rest ^ CS) A\{lab rest). Indeed, this transformation results in equal polar angles 
Ohi-h-z = arctan(gT/g). The inverse of this transformation takes vectors from the CS frame 
to the lab frame is: 



Ki'^iCS lab) = (A^^ilab CS))' 



\ 



g 

-— COS0y 
Qt ■ A. 

— — sm (pv 
^ 3 

Q 



Mt 

Mt 



V 



QMt 
cos (f) 

sin 0v 
QMt 





- sin (pv 

COS0V 





q' \ 



Mt 








Qrq' 
QMt J 



The kinematics of the leptons from the decay of the vector boson can be described by 
the polar angle 6 and the azimuthal angle 0, defined in the CoUins-Soper frame. The above 
transformation formulae lead to the four-momentum of the decay product fermion (and 
anti-fermion) in the lab frame as 
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where 




= {Mt coshy, Qt cos (pv, Qrsm0y, Mrsmhy), 
= ^ I q+n" + q^nf" - 

= €'"""^^Z^Xf3. (Al) 

Here, g± = j=^{q'±q'), y = ^Hq+ZQ-), n'^ = 73(1,0,0,1), = ^(1,0,0,-1) and the 
totally anti-symmetric tensor is defined as e^^"^"^ = — 1. 

APPENDIX B: 0{as) RESULTS 

To correctly extract the distributions of the leptons, we have to calculate the production 
and the decay of a polarized vector boson. The 0{as) QCD corrections to the production 



and decay of a polarized vector boson can be found in the literature in which both 
the symmetric and the anti-symmetric parts of the hadronic tensor were calculated. Such a 
calculation was, as usual, carried out in general number {D) of space-time dimensions, and 
dimensional regularization scheme was used to regulate infrared (IR) divergences because 
it preserves the gauge and the Lorentz invariances. Since the ant i- symmetric part of the 
hadronic tensor contains traces with an odd number of 75 's, one has to choose a definition 



(prescription) of 75 in D dimensions. It was shown in a series of papers |T6| that in D ^ A 



dimension, the consistent 75 prescription to use is the t'Hooft-Veltman prescription. Since 



in Ref. a different prescription was used, we give below the results of our calculation in 



the t'Hooft-Veltman 75 prescription. 

For calculating the virtual corrections, we follow the argument of Ref. and impose 
the chiral invariance relation, which is necessary to eliminate ultraviolet anomalies of the 
one loop axial vector current when calculating the structure function. Applying this relation 
for the virtual corrections we obtain the same result as that in Refs. and The final 
result of the virtual corrections gives 

as /47r/i2y 1 f 2 ^ , 2 o\ i ^^ |2 
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where e = (4— D)/2, /i is the t'Hooft mass scale, and Cp = 4/3 in QCD. The four dimensional 
Born level amplitude is 

\J^>-Born\ 



i9l + glM + fR)^o + (gl - - fl)As\ , (B2) 

where we have used the LEP prescription for the vector boson resonance with mass My and 
width Fy. The angular functions are Cq = 1 + cos^ 6 and ^3 = 2 cos^^. The initial state spin 
average (1/4), and color average (1/9) factors are not yet included in Eq. (P2|). 

When calculating the real emission diagrams, we use the same (t'Hooft-Veltman) 75 
prescription. It is customary to organize the ^(ct^) corrections by separating the lepton 
degrees of freedom from the hadronic ones, so that 



dQ^dydQld cos 6d^ J ^^^^^ {2nyS {Q^ - M^f + Q^VH 



with zi = Xi/^i and Z2 = X2/^2- The dependence on the lepton kinematics is carried by the 
angular functions 

£0 = 1 + cos^ 9, Aq = -{1 — 3 cos^ 9), Ai = sin 29 coscp, A2 = - sin^ 9 cos2<^, 
^3 = 2 cos 6*, Ai = sin 6* cos 0. 

In the above differential cross section, i = —1, ...,4 with A-i = Co', and = {gj^ + 5'^)(/l + 
/^) for i = —1, 0, 1, 2; ^* = {gl — gji){fl — /^) for i = 3,4. The parton level hehcity cross 
sections are summed for the parton indices a, 6 in the following fashion 

E ^o-b '^ah = E [f-qq Xjq + ^qq + ^qG TqG + ^qG TqC + ^Gq ^Cq + ^Gq '^Gq) ■ 

a,h q=u,d,s,c,b 

The partonic luminosity functions Cab are defined as 

Cab{^l,^2,Q'^) = fa/hiUl,Q^) fb/h2{^2,Q'^), 

where fa/hi is the parton probability density of parton a in hadron hi, etc. The squared 
matrix elements for the annihilation sub-process qq VG in the CS frame, including the e 
dependent terms, are as follows: 
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T^t = -AT^{u,t)-{t + ufe 



ut 



ut Ml 



For the Compton sub-process qG —>■ Vq, we obtain 

T-J = — (T4s,u)-{s + uye 
su ^ 

rSo = r^c = ((g^ - uf + (g^ + - (. + 



5 



XyG = — irr -2M(g - S) + 



= — ^ (2s(g2 - s) + T4s, u) - iQ' - t)ue) . 
^ su Mt ^ ' 

In the above equations, the Mandelstam variables: s = (k + l)"^, t = (k — q)'^, and u = {l — qY 
where k, I and q are the four momenta of the partons from hadrons hi, /12 and that of the 
vector boson, respectively. All other relevant parton level cross sections can be obtained 
from the above, and summarized by the following rules: 



-1,0,1,2 z = 3,4 



iTi '7~i n't n'i 

■^qG — ■'■qC -^qO ~ qG 

■^Gq — ■'■Gq ^Gq ~ ^Gq^ 



with the only exceptions that Tq^ = —T^g ^ ^) ^^"^ "^Gq ~ "^qG {u ^ t). These results 
are consistent with the regular pieces of the Y term given in Appendix E and with those in 
Ref. ra. 



In the above matrix elements, only the coefficients of Cq and ^3 are not suppressed by Qt 
or so they contribute to the singular pieces which are resummed in the CSS formalism. 
By definition we call a term singular if it diverges as Q^^ x [1 or ln(g^/g|,)] as Qt — >■ 0. 
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Using the t'Hooft-Veltman prescription of 75 we conclude that the singular pieces of the 
symmetric (Cq) and anti-symmetric (^3) parts are the same, and 

J™o 5is + t + u-Q^)= Jim^ T/^ 6is + t + u-Q') = s,^, 
lim Tag' 6{s + t + u-Q^)= lim 6is + t + u~Q^) = scg, 



where 



2 5(1 - 21) 5(1 - Z2) ( In ^2 



- ((1 - ^1) (5(1 - ^2) + (1 - ^2) 5(1 - ^1)) e] + O 



SGq 



y{4 + (1 - 4)) - ^2) - 5(1 - Z2)e\ +0\^-^ 



As Qt 0, only the Cq and ^3 helicity cross sections survive as expected, since the 0{ag) 
differential cross section contains only these angular functions [cf. Eq. (p^) ]. 

APPENDIX C: EXPANSION OF THE RESUMMATION FORMULA 

In this section we expand the resummation formula, as given in Eq. (|1|), up to 0{as), 
and calculate the Qt singular piece as well as the integral of the 0{as) corrections from to 
Pt- These are the ingredients, together with the regular pieces to be given in Appendix 0, 
needed to construct our NLO calculation. 

First we calculate the singular part at the 0{as)- By definition, this consist of terms 
which are at least as singular as Q^^x [1 or In {Qt/Q'^)]- We use the perturbative expansion 
of the A, B and C functions in the strong coupling constant as as: 



00 



B{as{-^),C,,C,) = Y,['^\ B'^^\C,,C2), (CI) 



n=l \ ^ 



00 



n=0 



The explicit expressions of the and coefficients are given in Appendix 0. 

After integrating over the lepton variables and the angle between b and Qt, and dropping 
the regular (Y) piece in Eq. (|I]), we obtain 
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da 



dQ^dydQl 



where we have substituted the resonance behavior by a fixed mass for simphcity, and defined 
ao as0 

"^o = ^T7^, fory = 7*, 



9Q2 



10„2 2 



Here a is the fine structure constant, sw (ciy) is the sine (cosine) of the weak mixing angle 
9wi Qj is the electric charge of the incoming quark in the units of the charge of the positron 
(e.g. Qup = 2/3, Qdown = —1/3, etc.), and Vjk is defined by Eq. (H). To evaluate the integral 
over rj = hQx-, we use the following property of the Bessel functions: 

/ dr]r] Jo{r]) F{r]) = - dr] r] Ji(r/) — - — , 
JO JO df] 

which holds for any function -^(77) satisfying [rj Ji^rj) F^t])]'^ = 0. Using the expansion of 
the Sudakov exponent S{b, Q, Ci, C2) = S^^^ (6, Q, Ci, C2) + 0(a|) with 

S^^\b,Q,C,,C.) = ^ [V)(COln^ (if ) (1^) 
and the evolution equation of the parton distribution functions 

we can calculate the derivatives of the Sudakov factor and the parton distributions with 
respect to rj: 



^For our numerical calculation (inside the ResBos Monte Carlo package), we have consistently 
used the on-shell scheme for all the electroweak parameters in the improved Born level formula for 
including large electroweak radiative corrections. In the V = case, they are the same as those 



used in studying the Z'^-pole physics at LEP |47|. 
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d 

dr] 
d 



,Q,CuC,) ^ -2 «5(Q') 
C|Q2 _2a5(Q2) 



ClQ^ 



+ and 



Note that as itself is expanded as 



27r 



27r 



/5o 



271 



In 



with Pq = {llNc — 2Nf)/Q, where A'',^ is the number of colors (3 in QCD) and Nf is the 
number of light quark flavors with masses less than Q. In the evolution equation of the 
parton distributions, 



and 



(C2) 



are the leading order DGLAP splitting kernels and ® denotes the convolution defined 
by 

(P^'la ® fa/,) (X, ^^') = £ I P^t fa/H 

and the double parton index a is running over all light quark flavors and the gluon. In 
Eq. (|C2|), the "+" prescription is defined as 

^ dz {G{z)). F{z) = dzG{z) [F{z)Q{z - x) - F{1)] , 
Jo 



where 



e{x) 



0, if X < 0, 

1, if X > 



is the unit step function and -^(2;) is an arbitrary function. 

After utilizing the Bessel function property and substituting the derivatives into the 
resummation formula above, the integral over 77 can be evaluated using 



dr]J,{7]) In'" 

\0n 



1, if m = 0, 

0, if m = 1,2 and 60 = 2e-^^, 



(C3) 



where je is the Euler constant. Fixing the renormalization constants Ci = C2&0 = C3 
60 = 2e~'^^, we obtain the singular piece up to 0{as) as 
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da 



dQ^dydQ^ 



T 



^(1) In f 



+ 



fj/hi{xi,Q^) /^/^^(X2,(5^) + j ^ k 



(C4) 



To derive the integral of the 0(05) corrections over Qt, we start again from the resum- 
mation formula [Eq. (|T])] and the expansion of the A, B and C functions [Eq. (|CJ1|)]. This 
time the evolution of parton distributions is expressed as 

fi/Ax, = f,i^(x, Q^) + ffl(x, t?) + 0(qD, with 



27r 



where summation over the partonic index a is implied. Substituting these expansions in the 
resummation formula Eq. (|I]) and integrating over both sides with respect to Q^. We use 
the integral formula, valid for an arbitrary function F{b): 



P2 



(27r)2 Jo 

together with Eq. ( p3|) to derive 



dQ^T / d'b e'^^^-'Fib) = — db PTJi{bPT)F{b), 



271 Jo 



p2 



X 



dQ 



do 



^ dQ^dydQl S 



V) 



TC 



2^ [Pl 



i?«ln 



PI 



(ysm, (Q 



2ti 



In 



p2 



[Pj^a ® fa/h,) {Xl, Q') h/h,{x2, Q^) 



vr 



+j ^ A; + / dQl Y{Qt, Q, Xl, X2) 



(C5) 



where xi = e^Q/\fS and X2 = e~^Q/\fS. Equations ( |C^) and ( |C5|) (together with the 
regular pieces, discussed in Appendix ^ are used to program the 0{as) results as discussed 
in the beginning of Section |T|. 
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APPENDIX D: A, B AND C FUNCTIONS 



For completeness, we give here the coefficients A, B and C utihzed in our numerical 
calculations. The coefficients in the Sudakov exponent are 0,^31 • 



A^^\C,) 



Cf 



67 
3 

'2 ~ 



TT 

2 In 

9 

TT 



Nc- 



Cf<CfI---- 3C(3) 



16 



36 



— + -C(3)) + 
48 2^^ ^ ' 



Nf 



-TT 



bo 



17 
12 



67 
18 



TT 

I HC2) 



c 



-Nf 



In 



+ 



where N j is the number of light quark flavors (m^ < Qy, e.g. N^ = 5 for or Z'^ 
production), Cp = tr(tata) is the second order Casimir of the quark representation (with ta 
being the SU(N(7) generators in the fundamental representation), [3q = (llA^c ~ 2Nf)/6 and 
({x) is the Riemann zeta function, and ({3) ~ 1.202. For QCD, Nc = 3 and Cp = 4/3. 
The C-^'' coefficients up to n = 1 are: 



Cjk{z,b,fi,^) 



6jk6{l - z), 
0, 

5jkCF\\{l 



+(5(1 



2 ^ 



In^ 



1 

Cf 
boC2 



-3/4 



7r^_ 23 
^ T ~ 16 



In ( ^\ P« 



bo 



where P'-]}_^{z) are the leading order DGLAP splitting kernels |^ given in Appendix C, and 
j and k represent quark or anti-quark flavors. 

The constants Ci, C2 and C3 = /i6 were introduced when solving the renormalization 
group equation for Wjk- Ci enters the lower limit /I = Ci/6 in the integral of the Sudakov 
exponent [cf. Eq. (^], and determines the onset of the non-perturbative physics. The 
renormalization constant C2, in the upper limit JI = C2Q of the Sudakov integral, specifies 
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the scale of the hard scattering process. The scale = C3/6 is the scale at which the C 
functions are evaluated. The canonical choice of these renormalization constants is Ci = 
C3 = 2e~'^'^ = 60 and C2 = C4 = 1 [0. We adopt these choices of the renormalization 
constants in the numerical results of this work, because they eliminate large constant factors 
inside the A, B and C functions. 

After fixing the renormalization constants to the canonical values, we obtain much simpler 
expressions of A'^^\ B^^\ A^"^^ and B^'^\ The first order coefficients in the Sudakov exponent 
become 

A«(Ci) = Cf, and B^^\Ci = bo,C2 = I) = -3Cf/2. 
The second order coefficients in the Sudakov exponent simplify to 

A^'\Ci = bo)=CF 



'67 vr2\ ^ 5^/ 



B^'\C,=bo,C2 = l) = C 



^ 4 
+CFNf 



^-3C(3)^ 

18 24 



^ \36 48 2^^ ' 



The Wilson coefficients Cj^ for the parity-conserving part of the resummed result are 
also greatly simplified under the canonical definition of the renormalization constants. Their 
explicit forms are 



b_o Ch 
b'C2 
b_o C, 

b' C2 



bo) = Sjk[^il-z) + ^5il-z){n^-8)^ and 



:z(l-z). 



As noted in Appendix 0, the same Wilson coefficient functions Cja also apply to the parity 
violating part which is multiplied by the angular function ^3 = 2cos^. 

APPENDIX E: REGULAR PIECES 



The Y piece in Eq. (|lD, which is the difference of the fixed order perturbative result and 
their singular part, is given by the expression 



Y{Qt, Q, x„ X2, 9, 0, Q) = ^ ^ £ 



TT 



xfa/hA^i, C,Q) R'^^'iQT, Q, zi, Z2, e, 0) h/h,{^2, C,Q) , (El) 
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where Zi = Xi/^i {i = 1,2). The regular functions R^^^ only contain contributions which 
are less singular than Q^'^x[l or ln(Q'^/Q'^)] as Qt — ^ 0. Their explicit expressions for 
hih2 — > V{-^ are given below. The scale for evaluating the regular pieces is C4Q. 

To minimize the contribution of large logarithmic terms from higher order corrections, we 
choose C4 = 1 when calculating the Y piece. 

We define the qq'V and the iii2V vertices, respectively, as 

ilt,[9L{l-l5)+9R{^ + l5)] and - Ts) + + 75)] • 

For example, for V = W~^,q = u, q' = d, £1 = u^, and £2 = e"*", the couplings g1 = fl = 
GfM^/\^ and 9%, = fn = 0, where Gp is the Fermi constant. Table | shows all the 
couplings for the general case. In Eq. (^), 
16 I 1 2 



R 



(1) 



l + gmi + WK + igl-gmi-mR'. 



where the coefficient functions i?"^ are given as foUow^^ 



QTn+{t,u) 



X 



( Q 



Q^KMt V^' QTMT'JZ+{t,u 



X 



[Ao + A2] 



Ai 



r: 



u s 



lQ_ f 2u{Q'-s) 



A, 



2Q2 



QtMj 



2s{Q^ - s) 



+ 1 



Aa 



with 



g2 (n+{t,u) 



Qt 



5{s + t + u-Q^) -2 5(1 - zi) 6(1 - Z2) 



-5(1 -^1) 



'l + z^ 

1-Z2, 



5(1 - Z2) 



l + zl 



01 



^'^Note that in Ref. [O] there were typos in and R2 
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and 

^G, = ^I-9t T^:±M + i + „ _ g2) _ M + (1 _ 

where TZ±{t,u) = {Q^ — t)^ ± (Q^ — u)"^. The Mandelstam variables s,t,u and the angular 
functions Co,Ai are defined in Appendix B. The Vjk coefficients are defined by Eq. (^). For 
a = j and b = G: = l^j^l^ where j and /c are light quark fiavors with opposite weak 

k 

isospin quantum numbers. Up to this order, there is no contribution from gluon-gluon initial 
state, i.e. Rgg ~ 0- '^^^ remaining coefficient functions with all possible combinations of 
the quark and gluon indices (for example R''^ , R'^^ or R^^ , etc.) are obtained by the same 
crossing rules summarized in Appendix 

Having both the singular and the regular pieces expanded up to 0{as), we can construct 
the NLO Monte Carlo calculation by first including the contribution from Eq. (|C5|) , with 
Pt = Qt'^^ for Qt < Qt^^- Second, for Qt > Qt^^, we include the 0{as) perturbative 
results, which is equal to the sum of the singular [Eq. (P^)] and the regular [Eq. pieces 
up to 0{as)- (Needless to say that the relevant angular functions for using Eqs. (|C^) and 
( PSD are Cq = 1 + cos^ 9 and A3 = 2 cos 9, cf. Eq. (p2|).) Hence, the NLO total rate is given 
by the sum of the contributions from both the Qt < and the Qt > Qt^^ regions. 
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